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by 

Ha rye 7 Brooks  and  P.  S.  Ham 
Harvard  "University 
Cambridge,  Massachusetts 

I. 

Introduction 

A number  of  papers1"1^  have  been  published  in  the  last 
twenty  years  concerning  the  determination  of  the  magnetic  sus- 
ceptibility of  a system  of  free  electrons  or  of  electrons 
confined  to  a box,  and  a variety  of  elegant  methods  have  been 
applied.  It  has  not  been  altogether  clear,  however,  to  what 
extent  the  walls  of  the  system  could  be  neglected  in  some  of 
these  calculations,  and  the  dependence  of  the  susceptibility  on 
the  form  of  the  wall  potential  has  not  been  adequately  investi- 
gated. Since  in  fact  the  walls  play  an  essential  role,  a deri- 
vation which  explicitly  considered  them  seemed  desirable.  The 
authors  believe  that  the  present  treatment  offers  a clear  picture 
. of  the  physical  situation,  some  worth-while  criticisms  of  results 
obtained  by  other  methods,  and  a few  useful  modifications  of  the 
- mathematical  methods  previously  used. 

It  is  well  known  that  in  classical  physics  the  positive 

moment  contributed  by  electrons  that  collide  with  the  walls  of 
the  system  exactly  cancels  the  negative  moment  of  electrons  far 
from  the  walls.  Van  Vleck^  and  Teller^  showed  that  a similar 
balancing  occurs  in  quantum  mechanics  (with  a non-zero  resultant) 
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and  gave  arguments  to  Justify  the  value  for  the  magnetic  suscepti- 
bility obtained  by  Landau^  from  the  free  energy.  However,  Teller's 
argument  was  for  an  infinite  plane  wall,  and  Van  Vleck's  required 
the  use  of  the  old  qnantum  mechanics,  and  it  is  not  clear  to  what 
extent1  either  Justifies  the  more  recent  calculations  by  Landau1 
and  others  of  the  oscillatory  De  Haas-Van  Alphan  terms.  We  shall 
first  give  an  improved  argument  for  a cylindrical  box  which  makes 
clear  this  balancing  of  diamagnetic  and  paramagnetic  states  and 
determines  with  the  WKB  approximation  the  nature  of  the  boundary 
states.  We  shall  then  give  a more  detailed  derivation  needed  to 
establish  corrections  to  the  terms  in  the  susceptibility  that 
dominate  at  high  field  strengths.  This  will  show  that  the  usual 
treatments1  give  the  De  Haas-Van  Alphen  terms  correctly  for 
"strong-  fields"  and  that  the  only  important  correction  to  the 

landau  steady  diamagnetism  is  a term  of  the  form  of  the  "surface 

7 8 

diamagnetism"  recently  reported  by  Osborne  and  Steele.  Unfor- 
tunately, the  magnitude  and  sign  of  both  this  surface  correction 
and  the  susceptibility  of  intermediate  and  small  systems  obtained 

•3  4. 

by  Dingle's  method^’  depend  rather  sensitively  on  the  assumptions 
made  about  the  form  of  the  wall  potential.  This  result  throws 
some  doubt  on  the  accuracy  of  the  WKB  method  in  such  calculations. 
The  following  derivation  is  valid  only  if  the  radius  R of  the  con- 
fining box  is  considerably  larger  than  the  classical  orbit  radius 
R = -%  l/2m£of  an  electron  moving  in  a plane  perpendicular  to  the 

C ,v  6X1 

magnetic  field  with  energy  equal  to  the  Fermi  energy  5 of  the 
system — the  "strong  field"  case.  Finally  we  shall  give  heretofore 
unpublished  results  for  the  model  proposed  by  Darwin10  to  show 
some  fundamental  differences  from  the  usual  electron-In-a-box  model. 

II. 

Methods  of  Calculating  the  Magnetic  Moment 

The  usual  procedure  has  ueeii  to  Calculate  the  magnetic 

moment  M from  the  Helmholtz  free  energy  F by  the  formula 
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(2.1) 


It  is,  however,  evident  that  if  the  electrons  are  independent  ex- 
cept in  so  far  as  they  obey  the  exclusion  principle  when  Fermi 
statistics  are  used,  then  the  moment  of  the  system  should  also 
be  given  by  the  sum  of  the  moments  of  the  individual  electron 
states  weighted  with  the  probability  that  the  state  is  occupied 
Thus  with  Fermi  statistics  we  should  have 


5,7 
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where  £ is  determined  by  the  condition 


N = 


? 


1+e 


(e1(H)-5)/kT 


(2.2) 


(2.3) 


N being  the  total  number  of  electrons  in  the  system,  and  the  sum 
being  taken  over  all  single  electron  states.  On  the  other  hand, 
the  free  energy  is  given  by 


F 


N 5 


(5-e.( H))/kT 
1 3, 


(2.4) 


and  since  by  (2.3)  = 0,  (2.1)  and  (2.2)  give  the  same  result 

provided  the  limits  of  summation  in  (2.4)  and  the  degeneracy  of 
states  with  a common  energy  (if  we  should  use  the  Index  i to  label 
a collection  of  states  with  the  same  energy  instead  of  a single 
state)  do  not  depend  on  H. 


However,  for  a system  confined  by  a cylindrical  box  of  cross 
section  A with  its  axis  parallel  to  the  uniform  magnetic  field, 
the  energies  of  states  unperturbed  by  the  walls  are 
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where  0 * is  the  component  of  the  propagation  vector 

along  the  axis,  and  the  degeneracy  of  each  such  (n,k„)  level  is 


8 «HA  a 

to  first  approximation0 * hc»  neglecting  spin  degeneracy. 

As  remarked  by  Van  Vleck^  and  Osborne,^  If  we  use  only  these 
energy  states  (2.5)  and  this  degeneracy  in  (2.2),  we  find  a 
large  negative  total  moment  of  magnitude  greater  than  N0,  whereas 
(2.4)aad  (2.1)  yield  the  usua,l  Landau  result (for(0H<< £) ), 


H - - 4afi.2J  (26&V3  H . _ fiaffl  • „ 

w2  K 9?;  a 2k  3 » 


(2.6) 


plus  the  periodic  De  Haas-Van  Alphen  terms.  Here  n(£)  is  the 
density  of  states  (including  spin  degeneracy)  in  the  energy 
scale  at  the  Fermi  level  £ , and  the  second  expression  shows  (2.6) 
to  be  much  less  in  magnitude  than  N0,  Our  first  problem  is.  then, 
to  show  that  the  difference  between  these  methods  js  removed  when 
we  include  the  positive  moments  of  the  boundary  states  in  our 
summation  (2.2). 


£he  Boundary  States  ^ the  WKB  Approximation 

We  consider  the  Schrodinger  equation  satisfied  by  the  wave 
function  of  a single  electron  state  in  a cylindrical  box  of 
length  L and  radius  R with  its  axis  parallel  to  the  uniform 
magnetic  field,  in  which  we  neglect  all  interaction  between 
electrons: 


- toLZ2  + $ + r ^(r  U + - E*  (3.1) 


for  r<  R,  and  the  boundary  condition  iKR) 
cylindrical  coordinates,  we  put 


On  separating  in 
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(3.3) 
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and  get  as  the  equation  for  f(r) 

' -jw  f=<r  Cr)  + + TM  - e’f(r). 

201  dr  dr  2mr2  8nc2 

We  now  make  the  substitution  r = ex,  f(r)  = g(x),  which  puts  (3.3) 
in  a form  suitable  for  the  WKB  approximation: 


afgisi  + 21  [E  ,ea.  fi|! . a!a!  e 

ax2  ii2  2“  8mc2 


4X 


]g(x)  * 0. 


(3.4) 


This  determines  e 1 in  terms  of  s and  n through  the  phase  integral 
condition  (after  transformation  back  to  r) 
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dr  = (n  + i)fc  tr, 


(3.5) 


where  rn 


and  r2  are  the  two  positive  zeros  of  the  quantity  under 


the  square  root,  provided  thar  r2  £ R.  This  condition  yields 

e1  = (2n  + | s | + 1)  3H,  (3*6) 

and  the  energies  are  given  by  (2.  5).  If,  however,  r^  < R <r2> 
then  the  upper  limit  of  integration  is  R,  and  it  is  evident  that 

1 n 

for  given  n and  s,  e'  will  be  increased  over  the  value  (3*6). 

We  note  that  the  quantum  numbers  n,  s,  and  k2  and  the  spin  orien- 
tation completely  specify  a single  state,  which  according  to  the 
exclusion  principle  can  be  occupied  by  no  more  than  one  electron. 

The  integral  in  (3.5)  with  R as  the  upper  limit  can  be 
1 ft 

evaluated  exactly,  but  the  result  is  a complicated  implicit 
equation  for  e 1 which  can  be  solved  only  with  great  difficulty. 
Portunately , when  s is  large,  the  effective  radial  "potential" 
in  (3*3), 


*2  2 2W2  2 
vCr)  = £-*3  + e_8-§- 
2mr<i 


8mc‘ 


(3.7) 


has  a sharp  minimum  at 
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and  to  excellent  accuracy  may  be  approximated  by  the  first  terms 
of  a power  series, 

v(r)  * pH | s ( + (r-r  )2.  (3.9) 

2mc 

Higher  terms  in  the  series  may  be  neglected  for  energies  of  the 
order  of  £ provided  R»  Rc,  the  condition  mentioned  in  the 
introduction.  For  the  present  we  consider  only  states  with  non- 
positive values  of  s,  so  that  we  have  from  (3.2),  (3*5),  and 
(3.9),  using  e * e1  - 0H|s|  « E - (ti2k2)/2m 


Setting  e = (2n+l)0H,  we  find  that  the  range  of  r for  which  the 
quantity  under  the  square  root  is  positive  is 


Ar  * 2r0 


(3.H) 


so  that  the  wave  functions  of  the  vast  majority  of  occupied  states 
are  highly  localized  radially  at  a distance  rQ  from  the  axis  of 
the  box,  provided  only  that 


eHR2  _ eHA 
2hc  he 


(3.12) 


the  number  of  states  with  common  spin  orientation  and  common  n and 
k_  and  rn  -^R,  is  very  much  greater  than  £/20H  ~ n+£.  This  is  in 
particular  true  for  boundary  states,  for  which  r ^ R.  The  quan- 
tity Ar  in  (3*11)  is  actually  the  diameter  -g  V2me  of  the  orbit  of 
a classical  electron  with  "transverse  energy"  e. 


We  determine  e as  a function  of  (rQ-R)  by  solving  (3.10)  with 
R replacing  r2  as  the  upper  limit  of  integration.  On  defining 

y-fc  . * - » (3-13) 


we  are  led  to  the  system  of  equations 
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4 . (2n+i)g(z)  , 

z2  cos  xz-  z(l-z2)* 


— «ffln)*VrT  = (ann)*f(*), 

(cos  ^z-zd-z2)*)*  - 


(3-14) 


where  z is  a parameter  with  the  range  -1  < z < 1.  These  have  been 
solved  numerically  for  n » 0*1,2,  and  the  results  are  plotted  in 
Fig.  I. 


The  moment  of  a sta^  of  given  n and  s is  given  by 


(3.15) 


For"bulk  states”  far  from  the  wall,  we  have  from  (3.10),  as  usual, 
e * (2n+l)0H  and  = -(2n+l)  0 , but  for  boundary  states  given 

by’ (3. 14)  it  is  readily  seen  that  the  first  term  in  (3.15)  is 
negligible  compared  with  the  second.  The  moment  is  therefore  posi- 
tive, since  (de  /dro)jj  > 0 and  from  (3.8) 


(3.16) 


Together  with  (3.13)  and  (3.14),  this  yields  for  boundary  states 


5u* 

0 


(2n+l)*(2  "tel  )*(l-z2)* 

i r 5-1 r-T  , C3.17) 

C0S_1z(C0S  AZ-Z(1-Z^)T)T 


which  is  plotted  in  Fig.  II  for  n = 0,1,2.  We  note  that  the 
moment  of  a single  boundary  state  is  proportional  to  the  square 
root  of  the  degeneracy  |s | Q given  in  (3.12).  As  we  shall  see  from 
(4.2),  (3.11),  and  (3.8),  the  number  of  occupied  states  is  propor- 
tional to  the  same  quantity,  so  that  the  total  moment  contributed 
by  boundary  electrons  is  proportional  to  |s|Q. 
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Approximate  Calculation  of  the  Total  Moment 


We  consider  first  the  system  discussed  by  Peierls,11  with  N 
electrons  in  a "two-dimensional"  system  with  common  value 
of  kz  and  common  spin  orientation,  at  the  absolute  zero  of  temper- 
ature. For  very  strong  magnetic  fields  the  degeneracy  •ft 
of  the  n = 0 level  is  larger  than  N,  so  that  all  electrons  are 
in  the  lowest  bulk  state,  and  the  moment  of  the  system  is  -N3*  As 
H is  lowered,  the  rQ  of  each  state  increases  by  (3.8)  until  some 
occupied  states  are  forced  up  into  the  "tail"  of  Fig.  I,  We  shall 
calculate  the  moment  of  this  tail  at  the  field  strength  at  which 
the  uppermost  filled  state  in  the  tail  is  Just  level  with  the 
n *=  1 hulk  states,  so  that  for  a further  decrease  in  H electrons 

must  "overflow"  from  the  n = 0 tail  to  n = 1 bulk  states,  if  the 

system  is  to  remain  in  thermal  equilibrium.  Then  from  (3.15). 
d 6 

neglecting  we  have 


where  p(rQ)  is  the  density  of  n = 0 states  in  terms  of  the  parameter 
rQ  and,  since  n and  s completely  specify  a state,  is  given  from 
(3;8)  by 

P(r0)  = (4.2) 


From  this  and  (3.16),  together  with  the  observation  that  for  a 
large  system  the  spread  in  |s|  of  the  states  in  the  occupied  part 
of  the  tail  is  negligible  compared  with  the  value  of  fs|Q  itself, 
we  see  that  we  can  take  |s|q  from  under  the  integral,  getting 


“tail 
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FIG.  I ENERGY  OF  THE  BOUNDARY  STATES. 
THE  QUANTITIES  x AND  y ARE 
DEFINED  IN  (3.13) 
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Since  the  filled  part  of  the  tail  stretches  between  the  unper- 
turbed bulk  states  n * 0,1,  we  have  finally 


‘tail 


H 


?pH  2N0 


(4.4) 


Hence  on  adding  to  this  the  moment  of  all  the  n c 0 bulk  states, 
whose  number  is  not  appreciably  diminished  from  N by  the  loss  of 
a few  states  to  the  tail,  we  see  that  the  total  moment  is  N0  , 
as  given  by  Peierls.  Continuing  this  argument  as  the  higher  bulk 
levels  fill  up,  we  find  the  same  oscillation  between  -N0  and  N0 
found  by  Peierls.  We  note,  however,  that  the  change  from  -N0 
to  N0  is  very  rapid  but  not  discontinuous  and  in  fact  occurs 
over  a range  AH/H  aj  a/VT7T0. 

It  is  now  a simple  matter  to  consider  finite  temperature  and 
to  justify  (2.1)  and  (2.4)  from  (2.2).  The  moment  contributed  by 
the  bulk  states  is  as  Usual  (for  given  spin  orientation) 


bulk 


eHA 

he 


00  00 
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whereas  the  moment  contributed  by  electrons  in  the  tails  is 

dr  p(r 


‘tails 


• h 


£ : / 
n=o  * od  tail  n 


1+e 


(4.6) 

Again  using  (3*16)  and  (4.2),  taking  |s|0  out  of  the  integral  on 
assuming  its  spread  negligible  over  the  part  of  the  tail  for  which 
is  appreciable,  and  using  (3.12),  we  have 


l+eu 
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On  carrying  out  the  integration  over  the  tails,  which  extend  from 
the  energies  of  the  unperturbed  n levels  to  infinity,  we  get 


6nll+e 


[£-(2n+l)3H 


.(4.8) 


Combining  this  with  (4.5)  we  see  that  we  have  for  the  total 
moment 


1 - 23  ? tn  (1+#l 

k n=o  J-  oo 


A*  „ 

C5-(2n+l)pH-  )^4.9) 


where  we  treat  £ as  a constant  parameter  in  the  differentiation. 
This  is  exactly  the  formula  obtained  from  (2.1)  and  (2.4)  when  we 
form  P using  only  the  energies  of  the  bulk  states  and  the  degen- 

S*cy  (3«12),  provided  we  use  (2.3)  to  eliminate  terms  involving 
and  neglect  states  of  positive  s.  Hence  to  the  extent  that 
s has  negligible  spread  over  the  occupied  portions  of  the  tails, 
the  usual  formula  (2.1)  is  valid  despite  the  field-dependent 
degeneracy.  We  now  recognize  the  vital  role  played  by  the  boundary 
states  in  cancelling  a large  part  of  the  bulk  diamagnetism,  and  we 
can  no  longer  suppose  such  calculations  applicable  to  electrons 
that  are  completely  free. 


Had  we  used  Boltzmann  statistics  instead  of  Fermi  statistics, 
the  arguments  would  have  gone  through  similarly,  and,  in  fact,  if 
we  change  our  parameter  of  integration  in  the  expression  analogous 
to  (4.7)  from  rQ  to  s,  we  get  essentially  the  boundary  term  in 
equation  (68)  of  Van  Vleck’s  derivation.^ 
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Correction  Terns  in  the  Total  Moment 

Instead  of  trying  to  refine  the  above  argument  directly  in 
order  to  obtain  correction  terms  when  the  spread  in  s is  appreci- 
able, we  now  turn  to  a more  elegant  derivation  which  is  more 
easily  handled  but  which  shows  the  nature  of  the  balancing  of 
bulk  dlamagnetlsn  and  surface  paramagnetism  less  clearly  than 
the  above.  We -consider  the  system  of  Section  III  with  the 
cyllndrlcally  symmetric  radial  potential  generalized  to  an  arbi- 
trary V(r)  so  that  in  the  Schrodinger  equation  (3*1)  BY  1*  re- 
placed by  (S-V(r))Y*  The  separation  of  variables  used  in  (3*2) 
is  still  valid,  and  we  now  define  rQ  by  (3*8).  Temporarily  we 
restrict  attention  to  non-positive  values  of  s and  require  that 
the  "transverse  energy"  parameter  e « E - (fi^kf/2m)  vary  suffi- 
cicntlyslowly  as  a function  of  |s|  to  allow  us  to  replace  sums  over 
|s|  occurring  in  expressions  for  the  moment  and  free  energy  by 
integrals  over  fel.  This  is  equivalent  to  the  requirement  that 
our  fields  be  strong,  that  is,  that  R»RC„  We  require  also  that 
the  length  L of  the  cylinder  parallel  to  the  field  be  long  enough 
so  that  sums  over  the  allowed  values  of  kz  may  be  replaced  by 

integrals  over  k_  with  the  density  L/2rr.  We  have  then  for  the 

* 19 

total  moment  of  electrons  with  a common  spin  orientation 
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(5.1) 


Integrating  by  parts  in  |s|,  noticing  that  contributions  from  the 
limits  of  integration  vanish,  taking  the  derivative  with  respect 
to  H (regarding  ^ as  a constant  parameter)  from  under  the  inte- 
gral signs  on  the  reasonable  assumption  that  d^/BHBs  ■ d2e/3sdH, 
and  finally  changing  our  parameter  of  integration  from  |s|  to  rQ 
oy  (3»6),  we  nave 


-11- 
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(5.2) 


The  bracketed  quantity  Is,  of  course,  the  quantity  (F-Jft)  (still 
neglecting  states  with  s positive).  We  may  also  derive  (5.2)  from 
(5.1)  by  first  changing  our  parameter  of  integration  to  rQ,  using 
(3.15),  and  then  performing  the  necessary  integration  by  parts. 

We  shall  now  consider  the  evaluation  of  (F-N§)  in  the  form  appear- 
ing in  (5.2). 


We  first  set  E ■ e ♦ 


iA? 


change  from  k to  E as  a parameter 


2m  ’ z 

of  integration,  and  integrate  by  parts  in  E,  obtaining 


F_N&  . i J £ J T^(  dr0> 

° “°°  (5.3) 

5(n,H,E)  is  that 
We  may  replace  S(n,H,E) 


where  f„  is  the  Fermi  function 


value  of  rQ  for  which  e 


E. 


l+e 


o 

by  oo  if  we  consider  (B-e)y  to  be  identically  zero  whenever  e>E,  as 
we  shall  do  henceforth. 


We  now  specialize  V(r)  to  be  constant  over  our  system  except 
for  a rapid  rise  at  the  wall.  Then  only  for  states  in  the  tail  will 
e depend  on  rQ,  and  we  may  replace  zero  in  (5»3)  *s  the  lower  limit 
of  integration  over  rQ  by  (R-A  rQn) , the  value  of  rQ  at  which  the 
nth  tail  starts  to  rise.20’21  The  tail  rises  rapidly  in  energy 
so  that  r2  differs  only  slightly  from  R2  for  occupied  states  in 
the  tall,  and  we  may  advantageously  use  the  expansion 


r2  = R2+  2R(r0-R)  ♦ (rQ-R)2 . (5.4) 

When  we  use  the  first  term,  R2,  in  (5*3) > the  integral  over  rQ  is 
immediately  and  exacUy  integrable.  The  only  contribution  comes 
from  the  lower  limit  of  integration,  where  the  energy  levels  are 

those  of  the  free  electron,  (2n+l)3H.  This  is  the  value  for  the 

1 2 3 

free  energy  derived  by  Landau,  Sondheimer  and  Wilson,  and  Dingle. 
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We  see  that  the  evaluation  of  this  term  In  the  expansion  is  ex- 
act, independently  of  any  approximations  that  must  be  made  in 
determining  £ in  the  energy  tail. 

The  evaluation  of  the  contribution  to  P of  the  second  and 
third  terms  in  (5.4)  will  depend  on  our  approximations  in  deter- 
mining e » In  the  following  work  we  shall  assume  the  V(r)  of 
Section  III,  an  infinite  potential  Jump  at  the  wall,  and  use  the 
approximate  results  of  (3.13)  and  (3.14).  Prom  these  equations 
we  see  that  to  the  accuracy  of  our  use  of  the  WKB  method  and 
the  approximate  potential  (3.9),  the  n^  tail  starts  up  at 
Aron  *(5f(2n41))  , that  (r0-H)  ■ ^^§(2n+l)^  * fCx) , and  that 
e * (2n41)0H  g(x).  Hence  it  is  convenient  to  change  our  vari- 
able of  integration  in  (5«3)  from  rQ  to  z,  which  together  with 
(5.4)  gives  us 


P-H5 


R2^(21g(n4i))*f(2) 

(5.5) 


+ (n+*)(f(z))2j  20H(n+*)  (E-23H(n+*)g(z))*. 

We  shall  now  consider  separately  each  of  the  three  terms,  which 
we  call  (P-N£)0 'vj  R2,  (P-N^^R,  and  (P-H$)2  l,  and  in  each 
we  shall  reverse  the  order  of  integration  over  z and  summation 
over  n.  We  see  that  we  have  to  consider  sums  of  the  form 


n=o 


(n+*)XOS-20H  (n+J-)g(z))y  , 


(5.6) 


where,  as  already  remarked,  we  consider  the  second  factor  identi- 
cally zero  when  2PH(n+£)g(z)  > E.  We  now  use  the  Poisson  sum 
formula^ 


CD  Qp  m 

2Z  «»♦*>  ■ 2l  (-1)r  / 


e 


2TTinr 


f(n)dn, 


(5.7) 
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according  to  which  (5*6)  equals 


y ' (-l)r(g)x+1  E7  j e2ni^ra-ux(l-u)y’du, 
r»-  ao 

where  we  have  put  £ - 2pfiiCz7  and  u * ^ I * 


(5.8) 


We  must  now  examine  the  function 


Wy) 


f elvu  nx(l-' 

Jo 


■u)7  du, 


(5.9) 


in  which  we  may  expand  the  exponential  in  its  Taylor  series  and 
use  the  beta-function  formula, 


(l-u)ndu  - r'<r?fJ+n+2?+li  > 


(5.10) 


to  obtain  the  expansion 


v - V7  iivii  r <jzm  . 

Vxy  ' / \ 3.'  r(x+y+s^2) 

s=o 


(5.11) 


Comparing  this  with  the  expansion22  of  the  confluent  hypergeometric 
function  m(z),  we  that 

V^cv)  - (iv)-i-“ 


(5.12) 


where  k = (y-x)/2,  and  m = (x+y+l)/2.  A useful  property  is  that 


dTxy(v) 

dv 


1TI+I,y(v)  • 


(5.13) 


The  asymptotic  expansion  of  V^yXv)  ^or  v > 0 is  derived  from  that 

of“k,m(z)s23 


y 9 -15“ 

V (v)^r(y+l)eiv  e 2 (y  ^ (v)-(y+1)Cl40(^) 

xy 

♦ r(x+D  cv)‘(xn)  (i ♦ o(^)). 


(5.14) 


The  ’expansion  for  v<  0 Is  obtained  from  this  by  replacing  v by 
its  absolute  value  and  talcing  the  complex  conjugate  of  the 
formula,  since  from  (5*9)*  ^^(-Ivl)  **  ^xy^v^’ 

Hence  (5*8)  becomes 


:2Bfig(*r)X+1  Txy(f>Hg& 


r=-  cd 


To  evaluate  (F-l£)0,  we  use  the  usual  procedure  on  (5.3)  in 
order  to  avoid  any  appearance  of  making  an  approximation,  replace 
ro  ky  Integra te  over  r0,  and  then  use  the  Poisson  sum.  The 
result  is 

CM®,  - $ <a)3/2  r « b?/2  z (-i)p  (?-i6) 

“ Jo  r»-  oo 

The  r * 0 term  is  evaluated  using  (5»9)  end  (5«10).  The  others 

TO 

may  be  changed  in  the  conventional  fashion  * with  the  formula 

Y ,Cv)  m 1 . 31  v ,(v)  * ± + -3L  - -2L  elv  V (-V).  (5.1?) 

v 2v  vo^v;  v 2v2  4v2  “*$° 


(5.1?) 


The  series  in  (5.l6)  from  the  first  term  in  (5.17)  vanishes  if  we 
pair  off  terms  with  r * +lrl  . In  the  series  from  the  second  term 
in  (5.17)  we  use  the  formula 


= . 4 , 

_2  12  » 


(5.18) 


s o tha  t 


(F-Bft, 


& *5/2 « 
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E*dE 


(5.19) 


- S T ^ £'  *w  «ff>; 

v/o  r*-  cd  r 

where  the  prime  on  the  sum  indicates  omission  of  the  r = 0 term. 

So  far  this  is  exact.  The  conventional  result1’2  is  obtained 
from  the  asymptotic  formula  (5.14)  (for  r > 0) 

T-*,0  <TP>~  ^ • (5.20) 

Higher  terms  in  (5*14)  lead  to  corrections  to  the  steady  and 
oscillatory  terms  of  (5*19)  which  contribute  negligibly  to  the 
corresponding  terms  of  the  magnetic  moment  provided  PH/£«1. 

In  evaluating  (F-H^  and  (F-l5)2  in  (5.5)  with  the  use 
of  (5.15)  f we  see  that  the  r = 0 term  in  the  sum  can  be  inte- 
grated exactly,  for  V^CO)  is  independent  of  z,  This  term  Van- 
ishes in  (F-H£)^,  but  not  in  (F-H£)2.  On  setting  z * sin9  for 
convenience  and  obtaining  g(z)  and  f(z)  from  (3.13)  and  (3.14), 
we  have  finally 


- 3 T 5“ 

Jo 

oo  , tt/J 

• £ (-1)r  / 
r=-  CD  J-T T/ 


E-5  dE 


(5.21) 


cos  0 sin  9 d9  V 


1 ol\PHg ( s inO ) J ’ 


2 X05na2  (pH)2 


T ^E7/2dE  V7!-!)'/72 

Jo  3B  ^ J—n/2 


(5.  22) 


sin29cos29d0  • 


2§i  g(sinO' 
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We  now  seek  the  contribution  to  P of  states  with  positive 
values  of  s,  whose  unperturbed  energies  are  found  from  (3.2)  and 
(3.6)  to  be  cA^/aa  *(2(n*s)+l) 0H.  If  R > 2RC,  states  per- 
turbed appreciably  by  the  wall  have  energies  above  the  Penal 
level,  so  that  we  may  calculate  the  contribution  to  P with  com- 
plete neglect  of  the  walls 

t 9*  / . I (5-(2(n+s)+l)  0H — *-*)/kT 

- -kT  ■&  / dk.  > > tali'*.  2“ 


ft 


z z 

n*o  s*l 


dE  > n(E-<n-*|)2(B)3/2 


(5.23) 


■*  I 

On  setting  n = n * - ■£,  we  can  use  the  Poisson  sum  (5.7),  end 

after  some  rearr.  sent  using  (5*17)  we  obtain 

(f-h§)+  - 31s  iMjgZ  - iflB  J 6& 


71  Jo  r=^ob  L1 


6 OH)  2 1»2 
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(5.24) 


We  may  now  proceed  most  simply  by  using  a proof  given  by 
Osborne^  that  for  an  electron  gas  in  a potential  field  of  cylin- 
drical symmetry  the  volume  in  quantum  number  space  bounded  by  a 
surface  of  constant  energy  and  the  plane  n * - $ is  independent 
of  the  field  strength,  at  least  to  the  accuracy  of  the  WKB  approxi- 
mation,24 and  that  consequently  if  we  replace  sums  over  quantum 
numbers  by  the  appropriate  integrals  the  resulting  approximation 
to  the  free  energy  is  independent  of  H.  In  the  present  discus- 
sion we  have  already  integrated  over  kz  and  replaced  y [ in  (5.1) 
or  ( 5*  2)  bv  9*  and  in  the  Pois  son  sum  formulo  the  r — 0 

Jen 


l^gKB 
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term  arises  when  f(n*»iO  is  replaced  b Y f(n+$-)dn.  We 

have  moreover  determined  onr  energy  eigenvalues  with  the  JKB 
method.  Consequently  if  in  (5*23)  we  replace  by  J ^ ds 

in  order  to  make  the  s integration  cover  the  range  (-cd,  ao) , 
we  should  expect  that  the  sum  of  all  the  r * 0 terms  should  be 
independent  of  H.  This  term  in  (F-W^)0  is  indeed  field-independent, 
and  this  term  in  (F-Hfc^  vanishes.  The  term  arising  from  the  men- 
tioned modification  of  (5.23)  is  the  first  term  in  the  first 
integral  in  (5.24);  which  combines  with  the  r * 0 term  in  (F-n£)2, 
the  first  in  (5*22).  The  sum  of  these  two  is  not  zero,  but  we 
recall  further  that  we  introduced  an  approximation  in  using  (3.9) 
for  the  effective  potential  and  found  accordingly  that  e was  a 


function  of  (rQ-R)  only.  The  exact  WKB  calculation  yields  a 
series  of  the  form 


«l(r0-B) 

,(r„-R)  + -J-B2 + • 


(5.25) 


Using  only  the  first  term  of  this  series  introduced  no  error  in 
(F-!?£)0,  but  use  of  the  second  term  of  (5.25)  in  (F-N£)^  yields 
terms  of  the  form  of  (F-Nj;)2.  The  r * 0 term  of  this  correction 
does  in  fact  cancel  the  r - 0 terms  of  (F-N£)2  and  (F-N£)+.  The 
derivation  of  the  second  term  in  (5*25)  is  outlined  in  Appendix  B 
for  the  purpose  of  obtaining  important  corrections  to  (F-N£)2» 
Osborne’s  theorem  spares  us  further  calculations  on  higher  terms: 
we  see  that  no  r = 0 term  will  contribute  to  the  magnetic  moment. 

Actually  we  obtained  (5.24)  from  (5.23)  by  direct  use  of  the 
Poisson  sum  without  changing  to  an  integral  over  s,  and  conse- 
quently we  obtained  the  second  term  in  the  first  integral  in 
(5.24).  This  is  got  rid  of  when  we  note  that  for  consistency  we 
should  have  used  the  Poisson  sum  in  (5«1)  instead  of  replacing 
\ 'f(s)  by  P f(s)  ds.  This  leads  to  replacing  the  sum  by 

*■8*  Xod  r*  f (s)ds , and  the  extra/*'  ds  f(s)  turns  out 

J-  QD  JO 

to  be  just  what  is  needed  to  cancel  the  second  term  in  (5.24). 
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The  contribution  to  the  magnetic  moment  of  (F  - and 

(F  - *£).,.  are  immediately  obtained  by  straightforward  differentia- 
tion with  respect  to  explicit  dependence  on  H,  with  use  of  (5*13)* 
From  the  asymptotic  expansion  (5*14)  we  find  that  the  dominant  steady 
and  oscillatory  terms  are 

9Paf, 


£ <|§)3/2  UOH)  f 3b2  b*  dE 


POH)* y** 


sin  (jj|®  - $),  (5.26) 
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(5.27) 
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To  obtain  the  contributions  of  (F  - n£)^  and  (F  - w©  must 

use  steepest  descent  procedures  to  evaluate  the  integrals.  This  is 
done  in  Appendix  A§  the  dominant  terms  of  the  result  are 
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where  £(5/3)  is  the  Riemann  zeta-function.  The  steady  term  of 
(5.29)  has  been  corrected  to  include  the  contribution  to  of 
the  second  term  in  (5.25),  the  derivation  of  which  is  sketched 
in  Appendix  B.  The  coefficient  of  the  oscillatory  term  of  (5.29) 
is  uncorrected,  for  this  term  is  too  small  to  be  of  any  real 
interest.  The  steady  term  of  M2  alone  is  3/2  times  the  steady 
term  of  (5.29). 

Tb  contrast  the  magnitude  of  the  various  contributions  to 
the  moment,  we  have  constructed  Table  I,  in  which  we  record  the 
form  of  the  dominant  steady  term  and  that  of  the  amplitude  of 
the  dominant  oscillatory  term  of  MQ,  11^,  M2+  Hls,  and  M^,  at 
the  absolute  zero  of  temperature,  where  df0/dB  c -8(E-§).  The 
correct  numerical  coefficients  are  given  only  for  the  steady 
terms  in  1IQ,  and  M2+  Hl8.  We  also  give  the  condition  on 
the  radius  R of  the  cylinder  and  the  field  strength  (Rc=  ^ V2m§) 
such  that  the  steady  terms  in  M^,  M2+  M^s,  and  M+  are  very  much 
smaller  than  the  steady  term  of  M0,  the  ordinary  Landau  diamag- 
netism. The  condition  given  for  the  oscillatory  terms  is  that 
required  to  make  them  small  compared  to  the  UQ  oscillations. 

Since  the  validity  of  our  analysis  is.  as  previously  remarked, 
restricted  to  "strong"  fields,  H»R  . we  see  that  the  only 

per  c 

significant  corrections  ' to  HQ  are  the  steady  term  of  M, , 
the  so-called  "surface  diamagnetism"  found  by  Osborne^’2  and 

Q 

Steele,  and  the  smaller  steady  term  of  “2+  “is-  The  numerical 
coefficient  and  temperature  behavior  of  both  are  given  in  (5*28) 
and  (5*29)  if  the  integrals  are  evaluated.  Of  course,  we  have 
given  only  the  dominant  terms  in  the  expansion  of  1£  asymptotic 
in  (jj/JJH).  Therefore  for  exceedingly  high  fields  such  that 
this  quantity  is  of  order  unity  (it  can  never  be  much  less  than 
unity  except  for  very  high  temperatures),  the  value  of  U will 
be  considerably  changed.  We  should  remark  that  all  of  our  results 
are  for  a system  of  electrons  with  common  spin  orientation.  They 
should  be  multiplied  by  two  for  a real  electron  gas  if  we  continue 
to  disregard  the  magnetij  moment  of  the  spin.  This  may  be  included 
by  the  method  of  Dingle.  We  may  thus  conclude  that  the  ordinary 


Table  I 


Dominant  Steady  Term 

Dominant  Oscillatory  Term 
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Landau  typo  of  calculation  using  the  free  energy  and  neglecting  the 
boundary  states  yields  the  correct  value  of  the  moment  provided 
R»  (5/fe  H)ly^Rc,  and  we  see  that  the  "surface  diamagnetism"  correc- 
tions may  be  obtained  in  a straightforward  manner  as  well  as  with 
the  elegant  number  theory  methods  of  Osborne  and  Steele. 


VI. 

Corrections  £&£  Sail  UalfijffiaH 


We  now  seek  the  effect  on  the  moment  of  the  system  of  replacing 
the  potential  used  for  detailed  calculations  in  Sections  III,  IV, 
and  V,  an  infinite  Jump  at  radius  B,  by  a potential  that  gives  the 
wall  region  seen  by  the  electrons  a finite  thickness.  Vs  assume 
the  potential 


V(r)  * a(r  - R)2  for  r > R 


s o 


for  r < R 


(6.1) 


and  in  order  to  simplify  our  work  require  that  when  aCr^  - R)2  = £ 
the  "thickness"  of  the  wall,  (r^  - R) , shall  be  very  much  less  than 
the  orbit  radius  Rc.  This  leads  to  the  requirement 


2tt2 


a» 


2mc 


2 * 


(6.2) 


”8  „7 

In  metals  the  wall  region  is  probably  of  the  order  of  10“  to  10"'  cm, 
and  Rc  is  of  order  10“^  or  10 for  fields  of  several  thousand  gauss, 
so  that  the  requirement  (6.2)  is  reasonable. 

We  must  now  include  the  potential  (6.1)  in  the  phase  integral 

(3.10): 

R T. 
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[2m(e - (r-r  )2)]*  dr  + [2m(e.  (r-r  )2-a(r-R)2)  ]*  dr 
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(6.3) 


Changing  variables  with  (3«13)  and  setting  u = (rQ  - R)C^)^  9 
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_a2L 


2m  3 H 
involving  y 


tJ“2  , and  noticing  that  because  of  (6.2)  y»l,  we  obtain  on 
integrating  and  keeping  only  the  largest  of  the  terms 


t.2n+3> 
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x2  4 

where  as  usual  y = to  order  l/yT,  and  -1  < z < 1< 


We  now  recall  that  our  evaluation  of  (P  - n£)q  in  Section  V did 


not  depend  on  the  particular  shape  of  the  energy  curve  plotted  against 
rQ,  so  that  a potential  of  the  present  sort  in  no  way  changes  the 
previously  obtained  value  of  (F  - n£)q  or  Mq.  The  only  corrections 
come  in  through  (P  - N&)^  and  (P  - N&Jg.  We  shall  consider  only  the 
former,  the  larger  one.  We  must  then  use  (6.4)  and  (6.5)  in  the 
integral  (B4)  of  Appendix  B>  We  use  the  same  device  outlined  there 
to  put  the  integral  in  standard  form  and  use  the  methods  of  Appendix 
A,  with  6W  = cos^O  , to  evaluate  it.  The  r = 0 term  does  not 
vanish  but  is  independent  of  H,  as  we  expect  from  Osborne’s  theorem. 
The  oscillatory  part  of  this  correction  to  the  moment  is  entirely 
negligible  for  the  conditions  we  have  assumed.  In  addition  to  the 
usual  steady  "surface  diamagnetism”  correction  (5»28),  however,  we 
have  the  diamagnetic  steady  contribution 
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This  turns  out  to  be  small  compared  to  the  usual  "surface  dia- 
magnetism" (5.28)  only  if 
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This  is  a stronger  condition  than  (6*2),  but  one  that  is  satisfied 
in  practice  with  a factor  of  10^  or  so  to  spare,  using  the  estimates 
for  Bc  and  wall  thickness  cited  above.  Hence  Mlw  is  negligible 
compared  to  We  may  infer  that  this  conclusion  is  true  for  any 

other  form  of  wall  potential  confined  to  an  appropriately  thin 
layer. 

The  above  considerations  make  it  appear  that  the  susceptibility 
is  at  least  for  "strong"  fields  relatively  independent  of  the  details 

t 

of  the  wall  potential.  This  is  true  for  Mo,  since  the  calculation  , 

of  this  term  was  exact,  but  it  is  not  true  for  M^  for  the  following 

reason.  In  applying  the  WEB  phase  integral  (3.5)  to  determine  the 

energy  eigenvalues  we  should,  as  remarked  in  a footnote,  ^ replace 

(n+£)  by  (n+3/4)  when  we  assume  an  infinite  potential  wall  (or 

equivalently  the  boundary  condition  t(R)  = 0)and  when  the  upper  limit 

of  integration  in  (3.5)  is  B.  Such  a replacement  is  dictated  by 

WEB  theory  in  the  requirement  that  the  phase  of  the  approximate  wave 

function  in  the  classically  allowed  region  between  the  turning  points 

be  such  that  the  wave  function  vanishes  at  the  turning  point  instead 

of  connecting  to  a damped  exponential  in  the  classically  forbidden 

on 

region  beyond  the  turning  point,  the  more  usual  condition.  Similarly, 
if  both  turning  points  occur  at  infinite  walls,  we  must  use  (n+1) 
in  (3.5)  with  n = 0,  1,  2,  ...  as  usual.  This  latter  form  of  (3.5) 
in  fact  yields  the  correct  eigenvalues  for  a one -dimensional  electron  i 
in  the  potential  V(x)  = 0 for  |x|<a,  V(x)  = oofor  |x|>a.  The 
(n+3/4)  form  of  (3.5)  yields  the  exact  eigenvalues  for  a similar  j 

problem  with  V(x)  = oo  for  x<0,  V(x)  = kx2  for  x>0.  Consequently 
if  we  assume  an  infinite  potential  Jump  at  the  surface  of  the  con-  J 

tainer  in  the  present  problem,  we  must  use  (n+3/4)  in  (3° 5)  when  ; 

the  upper  turning  point  is  at  R.  ! 

1 

At  first  sight  one  would  not  expect  this  change  to  make  much 
difference,  since  n is  quite  large  for  most  of  our  levels.  However, 
it  turns  out  that  the  change  drastically  alters  the  Poisson  sum 
formula  (5-7)  and  in  this  fashion  alters  the  surface  correction  terms 
and  Mg  in  such  a way  as  to  produce  a "surface  paramagnetism" 


j 

4 


! 


TH169 


-24- 


instead  of  the  "surface  diamagnetism"  found  earlier.  The  para- 
magnetic result  agrees  with  Dingle's  calculations $A  in  fact  our 
result  for  and  X2  using  (n+3 A)  agrees  numerically  with  the 
terms  of  corresponding  form  and  magnitude  in  Dingle's  result 
(when  ours  are  multiplied  by  2 to  include  spin  degeneracy).  This 
is  as  it  should  be,  for  although  Dingle  begins  his  derivation  by 
considering  the  location  of  the  zeros  of  the  approximate  WEB 
wave  functions,  his  fundamental  equations  are  more  easily  derived 
by  using  (3.5)  with  (n+3A)  to  determine  the  number  of  eigenstates 
of  given  s and  kz  with  energies  below  any  given  E,  the  turning 
point  for  such  s,  k , and  E being  at  R. 

The  modified  form  of  the  Poisson  sum  formula  is  derived  in 
Appendix  C.  If  0 < a < 1, 


y e-3rira  F 

n«o  r=-  oo  •© 


(6.8). 


Using  this  in  place  of  (5. 7),  we  find  we  must  replace  the  (-l)r  in 
the  sums  of  (5«21)  and  (5.22)  by  e“^rr^Tat  Calculation  of  the 
resulting  numerical  coefficients  for  and  M2  + for  a number 
of  values  of  a between  1/2  and  3/4-  yields  Table  II,  in  which  the 
coefficients  replace  the  bracketed  numerical  coefficients  of  Table  I. 


Hence  the  magnitude  and  sign  of  the  surface  corrections  depend 
critically  on  the  choice  of  a.  The  same  is  true  for  the  value  of 
the  susceptibility  at  weak  fields,  R*£<RC,  calculated  by  Dingle.4 


Table  U 


a 

Coefficient 

Coefficient 

.75 

+ .00763 

- .00139 

.70 

- .00047 

+ .00009 

.65 

- .00770 

4 .00140 

o 

NO 

o 

- .01378 

+ .00251 

. **?  1 

• + + 1 

- ,01896 

+ .001 4. 

- ■ - • + 

.50 

- .02293 

+ .004173 
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As  may  be  seen  from  his  derivation  the  coefficient  of  this  term 
depends  on  the  form  of  the  Poisson  sum  exactly  as  do  the  coeffi- 
cients of  Table  II  and  would  yield  a paramagnetism  roughly  three 
times  the  magnitude  of  the  reported  diamagnetism  if  a e 1/2  were 
used.  All  this  seems  contrary  to  what  our  physical  intuition 
would  lead  us  to  expect,  and  we  should  be  Justified  in  maintain- 
ing a healthy  skepticism  toward  the  accuracy,  or  even  the  reality, 
of  surface  corrections  and  weak  field  terms  calculated  in  this 
fashion  from  the  WKB  approximation.  Such  skepticism  is  especi- 
ally Justified  if  we  consider  the' present  work  to  be  an  accurate 
approximation  to  the  state  of  affairs  in  an  actual  metal,  where 
the  surface  is  not  sharp. 

However,  if  we  assume  the  reality  of  the  surface  correction 
and  weak  field  susceptibility  and  the  accuracy  of  the  VKB  method 
in  calculating  them  when  we  assume  a general  wall  potential  more 
closely  approximating  the  wall  potential  seen  by  an  electron  in 
a real  metal,  we  must  decide  what  value  of  a to  choose.  Cer- 
tainly the  wall  is  not  infinitely  sharp,  so  that  we  cannot 
& priori  use  a ■ 3/4,  but  it  may  be  sufficiently  sharp  so  that 
a * 1/2  is  not  sufficiently  accurate  for  present  purposes.  To 
attempt  to  establish  a criterion  for  choosing  the  value  of  a we 
may  consider  a one-dimensional  electron  bound  by  the  potential 


V(x) 


C co 

0 


x < -L 
-L  < x < 0 
x > 0 


(6.9) 


On  one  hand  we  shall  obtain  the  exact  solution  for  the  energy 
eigenvalues,  on  the  other  the  approximate  HKB  solution.  Compari- 
son of  the  two  will  allow  us  to  choose  the  best  value  for  a. 

The  exact  solution  for  the  wave  function  for  -L  < x < 0 and 
the  value  of  the  logarithmic  derivative  at  x = 0 are 


t(x)  * A sin  k (x  + L) 
■ k cot(kL) 


(6.10) 
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where  k * ^ (2mE)^.  Por  x > 0,  we  must  solve  the  equation 

E*  (6.11) 


subject  to  the  boundary  condition  that  $(x)  approaches  zero  as 
x-*+od.  Making  the  substitutions 


a> 


mX 

H2 


E 


+(*)'  uCfe) 
we  obtain  the  differential  equation 

,2s  rt'  dl  (§) 


* (n  + £)Hiw 
£ * ax 


(6.12) 


d|2 


2| 


JL 

dfc 


2n  B^(6) 


(6.13) 


This  is  the  differential  equation  of  the  Hermite  polynomials  when 
n is  integral,  but  we  must  here  consider  a more  general  n.  We 
generalize  from  the  well-known  formula  for  the  Hermite  polynomials 


(-1)"  e*‘ 


afe,n 


(-ire 


Sr"'/ 


— z 


(*-|) 


n+1  dz  i 


(6.U) 

where  the  contour  circles  the  point  z = % once  counterclockwise.  We 
find  that  for  arbitrary  n the  function 

.2 


2^2 

2frl 


f- 


(z-|) 


n+1 


dz 


(6.15) 


satisfies  (6.13)  provided  the  integrand  has  the  same  value  at 
the  beginning  and  end  of  the  contour.  This  condition  and  the 
boundary  condition  u^(^)  ► 0 as  £ — + cd  are  sat- 

' r, 

isfied  if  we  cut  the  z-plane  along  the  real  axis  from  z to 
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start  the  contour  at  +w  on  the  upper  side  of  the  cut,  loop 
once  counterclockwise,  and  return  to  +cd  below  the  cut* 

It  is  now  easy  to  evaluate  ^(0)  and-  by  substituting 


t * The  contour  D in  the  t- plane  p-o  corresponding  to 

C makes  use  of  two  Riemann  surfaces  in  looping  the  origin  twice 
but  it  may  be  deformed  into  two  loops  around  the  origin,  each 
stretching  out  to  ♦co  on  the  real  axis.  We  may  now  replace  this 
contour  by  C*,  which  has  the  form  C had  in  the  z-plane,  and  we 
now  use  a single  Riemann -surface  eut-froathe  origin  to  +oo 
along  the  real  axis  but  traverse  C1  twice,  the  phase-of  the 
Integrand -on  the  second  circuit  being  advanced  over  that  of  the 

I 

first-circuit -by  the . difference-  in  phase  of  the  Integrand  on  the 
two- sides  of  the -cut*-  The  value  of.  the  resulting  contour  integral 
is  now  obtained  immediately  from  Hankel’s  integral  for  the  gamma- 

pO 

function. 


V"> 


u<°> 


-nin  m 
e cos  -s? 

r(l  ♦ n/2) 

2e“n^nsin 
P(i  ♦ n/2) 


(6.16) 


We  may  therefore  obtain  the  logarithmic  derivative  of  the  wave 
function  at  the  origin 
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Ha  (0) 

(oT 
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dl 

ds 


tan  in  SSlVJ 

2 rt*  ♦ h 


(6.17) 


•2  2‘ 

Equating  this  to  (6.10),  we  obtain  the  equation  for  E whose  roots 
yield  the  eigenvalues? 

“ + 1) 


k cot  kL  * 2(^f)1/4  tan 


p<i 


ref  + \ ) 


(6.18) 


For  the  WKB  solution  we  use  the  phase  integral  condition 
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where  s is  a positive  integer  and  we  consider  E to  be  some  root  of 

(6.18)  and  a to  be -some  number .between  zero  and. one  which  makes 

(6.19)  exact.  Integrating  (6.19)  and  using  (6.12),  we  find 

(s+a)  tt=  kL  + (2n+l)  * . (6.20) 


We  now  use  (6.20)  and  (6.12)  to  find  that 

cot  kL  = cot  | (s+a)  tt-  WJ  = cot  | 
Substituting  this  in  (6.18),  we  obtain 

r(-£r.  * 4) 


rra  - 


.(6.21) 


c<*{na  ' ic}  = (^)i  r(j£  ~ I cot{^  *(6-22) 

1 4' 


We  now  assume  that  L is  large  enough  so  that  the  spacing  of  the 
eigenvalues  of  the  system  is  small  compared  to  tiu)  , so  that  we  may 
regard  E in  (6.20),  which  really  denotes  an  eigenvalue,  as  a variable 
with  a continuous  range.  We  then  solve  (6.21)  for  a as  a function  of 
E/fltoo  Values  at  representative  points  are  given  in  Table  III. 


mas  Ill 


E/fto 


a 


1 

1/2 
1/10 
1/100 
-*•  0 


-►3A 

.750 

.753 

.743 

.750 

.838 

.938 


We  see  that  provided  E/ha)>l/2.  a is  remarkably  close  to  3/4. 
We  may  therefore  expect  for  the  WKB  solution  to  the  problem  of  the 
electron  in  a box  in  a magnetic  field  that  a will  be  close  to  1/2 
provided  E/fua  >1/2,  since  in  this  instance  the  potential  is  slowly 
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i 


varying  at  the  other  turning  point  in  contrast  to  the  above  model. 
Approximating  the  wall  potential  in  an  actual  metal  as  in  (6.1)  by 
-J-XCr  - R'2,  and  defining  X in  terms  of  the  wall  thickness  (rx  - R) 
and  the  Fermi  energy  by  i X(rx  - R)2  * *e  find 


(rT-R) 
■.  «-8 


(^)* 


(6.23) 


Dslng^lO"12  erg.  (Pj  - B)~10*8c»,  and  the  free  electron  mass, 
we  find  sjyv  1/10.  Comparing  this  with  Table  III  we  see  that  this 
puts  us  in  the  critical  region  with  an  intermediate  value  of  a, 
especially  since  most  occupied  states  have  energies  below  £ . How- 
ever. we  see  from  this  that  £ m has. to  be  quite  small  before  we  get 
into  the  "sharp  wall  region"  where  a * 3 A ia  magnetic  problem, 
and  we  see  from  Table  II  that  in  the  intermediate  region  the  dominant 
surface  correction  is  diamagnetic  (and  Dingle's  weak  field 
susceptibility  is  paramagnetic) , although  not  so  large  as  we  found 
previously  using  a «*  1/2.  This  suggests  that  provided  the  periodic 
potential  doesn't  alter  the  situation  too  drastically  in  a real 
metal,  the  surface  correction  is  diamagnetic.  However,  skepticism 
of  the  accuracy  of  all  these  susceptibility  calculations  seems 
justified. 

Finally,  we  may  remark  that  the  theorem  proved  by  Osborne 

(Section  V)  requires  slight  modification  when  (n+£)  is  replaced  by 

(n+a)  in  the  calculation  of  the  energy  eigenvalues  by  the  WKB  method. 

As  may  be  seen  from  Osborne's  proof,  the  volume  in  quantum-number 

space  that  is  invariant  is  now  bounded  by  the  plane  n = -a  instead 

0f  n = Fortunately,  this  is  the  natural  generalisation,  since 

Jthe  modified  Poisson  sum  (6.8)  leads  us  to  replace  f(n+a)  by 
fXB  n=0 

f(n+a)dn  in  the  r = 0 term. 


▼II. 


Tflfi  Harmonic  Oscillator  Model  fif  Paffltln 

We  have  heretofore  required  that  V(r)  be  constant  except  for  a 


i 
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rapid -rise  at  the  wall*  That  the  results  can  be  quite  different 
when  the  boundary  is  extended  is  shown  by -the  following  results  on 
Darwin’s  models  These  are  of  no  importance  in  the  interpretation 
of  experimental  results,  as  the  model  is  rather  unrealistic,  but 
they  are  given  for  curiosity's  sake,  since  the  problem  seems  to  be 
the  only  one  that  can  be  solved  more  or  less  exactly  for  all  ranges 
of  field  strength,  as  well  as  to  show  that  contrary  to  intimations 
in  the  literature  Darwin's  model  does  not  give  the  same  results  as 
the  electron- in-a-box  model  when  Fermi  statistics  are  used.  Hence 
it  should  be  used  with  care  if  used  at  all  as  a guide  to  thought. 

Following  Darwin,^®  we  consider  a system  of  electrons  bound 
by  the  potential  of  the  two-dimensional  harmonic  oscillator: 

Y = mw2(x2+y2).  We  first  use  an  approximate  procedure  for  strong 

fields  to  show  some  interesting  features.  We  choose  to  such  that 
fttu  «0H  for  H > Hq.  Then  we  may  show  that  with  neglect  of  the 
spin  the  energies  of  individual  states  are  to  a sufficiently  accurate 
approximation  for  H > Ho 

2 2 *2*2  0 < n <oo  , 

en ,m',k  ’ <2n+l)ffl  + m* ^ + -g*  , co  , (7.1) 

where  n,m',kz  and  spin  orientation  form  a complete  set  of  quantum 
numbers . 

If  now  we  consider  a "two-dimensional"  system  of  electrons 
with  a common  value  of  kz,  it  proves  quite  easy  to  sum  the  moments 
of  the  lowest  N states  to  get  the  total  moment  of  the  system  at  ' 

absolute  zero.  This  is  sketched  in  Fig.  Ill  as  a function  of  the 

p ! 

parameter  y related  to  H by  i 

I i 

h2  . y2  # (7.2)  . I i 

2p 

This  is  qualitatively  different  from  the  results  for  the  electrons  in 
a box  in  that  the  system  never  becomes  paramagnetic  and  the  ampli- 

O * 

A J _ _ A>  J.  ...  . ...  _ _ TT  J.  UAn/v/.../Mt  //  1 

u uue  yjx  usu  ilia  oxun  ucv.ioaoco  as  it  ucvtca.coo  — uu>vt  , tv*  I ' - » j 

the  local  average  of  the  moment  is 


FI6.ni  MAGNETIC  MOMENT  OF  A TWO-DIMENSIONAL 
SYSTEM  OF  ELECTRONS  IN  OSCILLATOR 
POTENTIAL  AT  T=0°K.  SPIN  IS  NOT  IN- 
CLUDED. (SEE  EQUATION  (7.2)  OF  TEXT  FOR 
DEFINITION  OF  15 2 
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T=  & , (7.3) 

3 

i 

where  n(£)  represents  the  density  of  states  at  the  Fermi  level,,  In 

the  box  problem  the  average  moment  is  zero  for  the  two-dimensional 

model o We  can  add  the  average  moments  (7-3)  for  systems  of  electrons 

with  different  values  of  k and  get  the  result  in  the  same  form,  with 

z 

n(£)  the  total  density  of  states  at  the  Fermi  level  of  a three- 
dimensional  system.  This  is  now  formally  exactly  the  same  as  one  of 
the  forms  of  Landau* s result  for  free  electrons  in  three  dimensions 
given  in  (2.6).  Of  course,  n (£)  is  different  in  the  two  problems  and 
in  particular  depends  on  to  here.  The  harmonic  oscillator  model  also 
leads  to  low  temperature  oscillations  in  the  susceptibility  per 
electron  of  the  complete  system,  as  plotted  for  absolute  zero  in 
Fig.  IV.  Again  the  curve  differs  from  that  for  the  system  of  electrons 
in  a box,11’12  in  having  the  oscillations  smooth  and  much  smaller  in 
magnitude  than  the  Landau  average  value.  In  Fig.  V the  susceptibility 
is  plotted  when  the  contribution  of  the  spin  to  the  moment  is  in- 
clude dc 


The  curves  of  Figs.  IV  and  V were  plotted  from  numerical  calcu- 
lations at  absolute  zero.  We  now  turn  to  a new  derivation  for  finite 
temperatures  and  arbitrary  field  strengths,  which  will  be  a slight 
modification  of  the  method  recently  introduced  by  Sondheimer  and 
Wilson  (S-W).2 

Following  S-W,  we  seek  first  to  evaluate  the  classical  partition 
function  Z(Y),  where  y c 1/kT* 

TT — 7 “Y®4 

z(y)  =2.  e 1 • (7o4) 

i* 

The  Schrodinger  equation  for  the  Darwin  model  is  (3.1)  with  E 
replaced  by  (E  - -Jm  o>  2(x2+y2) ) . We  let  our  system  have  length  L 
in  the  direction  parallel  to  the  field.  The  energy  eigenvalues  are 

E _ . * .(2n  + |s|+l)(32H2+fi2o>2)*  + s3H  + — , (7.5) 

n?S9KZ 
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n = 0,1,2,.. ..  j s 0 , ^1 j +2 ^ « • • • • 

The  approximate  values  (7.1)  may  be  obtained  from  (7*5)  when 
Ao>  «3H. 


We  assume  L large  enough  so  that  we  may  replace  the  sum  over 
k2  by  an  integral.  Then 


Z(r)  = 


e 


-rL2n+|  s|  +1)  T +s  pH]  £_ 

2n 
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A^k. 


_v  — 

’ 2m 

e dk. 


L(^=|)* 

Th2 


(7.6) 


e 


-rC2n+|s|  +l)?+s3H] 


n*o  s*-oo 


where  T = ( 32H2+A2go2)^  . 

The  sums  over  n and  s are  easily  done,  and  the  result  is 


Z(Y)  = ^ (-22B)i 1 

4 yn  sinhC^  (T-$h)  3sinh[-|cr+3H)  ] 


(7.7) 


It  is  an  interesting  verification  of  Dingle  and  Osborne's  theorem 
(Section  V and  Reference  24)  that  if  in  (7*6)  we  replace  sums  over 
n and  s by  integrals,  the  integration  limits  for  s being  (-oo,  oo), 
and  those  for  n (-  i « oo),  then  Z(y)  turns  out  to  be  independent 
of  H.  If  the  n lower  limit  nQ  is  anything  other  than  - *|,  however, 
Z(y)  depends  on  H through  the  factor  e“YT(2n0fl)^  vaiue  „ I 

occurs  automatically  in  the  Poisson  sum  formula  and  is  required  in 
Osborne's  proof. 

We  now  use  directly  S-W's  conclusion  that  the  free  energy  for 
Fermi  statistics  is  given  by 

OD  ^ j - 

F * N£  = f 5S2  z(E)dE , (7.8) 
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z(E) 


C+i  00 

r jx  aai  dK 

4 ^ 

e-i  oo 


(7.9) 


the  path  of  Integration  being  to  the  right  of  all  singularities 
of  the  integrand,,  As  before,  our  results  are  for  electrons  with 
a common  spin  orientation,. 

The  integral  (7.9)  is  evaluated  in  a straightforward  manner 
by  the  method  of  S-W.  As  in  S-W,  the  poles  of  the  integrand  con- 
tribute the  oscillatory  part  of  z(E) , the  integral  around  the 
branch  point  at  the  origin  the  steady  part.  The  general  result 
is 


z(E) 


f.2EI2h£ 


a=l,2  n^ 


(“Dn 

(2nn) 


rr+(~i)aP  h)3/2 

s in  [S§-5- & - ( -1 )' apH) 2 3 


• sin  ( — ^rm^— — “ H) 

r+  <,-i)apH  * 


(7.10) 


»7/2  b2h2  k3/2 

+ A 1 J - -^T  (1  + 2 -*pp)  + ....  , 

1&2P(9/2)  12  fi2^2  T(5/2) 

where  higher  terms  in  the  asymptotic  expansion  of  the  steady  part 
are  small  provided  E 'pH  » 1,  E/ftu»  1, 

For  high  fields,  PH  » tlco,  we  can  most  simply  obtain  the 
dominant  terms  of  (7.10)  by  replacing  sinh  ^(T-Pl)  in  (7»7)  by 
•^CT-pH)  and  integrating  (7. 9)  anew.  Using  this  result  to  calcu- 
late the  moment,  as  in  S-W,  we  find  with  neglect  of  small 
tempeiature-deoendent  corrections  that  the  dominant  terms  are, 
in  terms  of  n(£) , the  density  of  states  at  the  Fermi  level, 


_ /TTiS  n, 
r cos(q-tj-  - t, 


sinh( 


(7.11) 


which  may  be  contrasted  with  the  usual  result^-  for  electrons  in 
a box,  with  the  appropriate  n(£) 


I 
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llL6i£  * „,(.v>2_/kTx,_§  >4^  C_11r  sln(m*-  4' 


sinh( 


(7.12) 


The  principal  differences  are  in  the  phase  and  amplitude  of  the 
oscillations.  Thus  the  amplitude  in  (7.11)  varies  as  (-&)£, 
whereas  that  of  (7.12)  varies  as  the  reciprocal  of  this  and  is- 
therefore  much  larger  than  the  Landau  steady  term  for  suffi- 
ciently low  temperature. 

^fields  low  enough  so  that  £J  «1,  but  large  enough  so 
that  We  obtaln  fro“  <7.10)  the  dominant  terms 


i-.fl&fi?. 


tvr  n(§)32(f>*  ##  • 

•^b  co3<^^‘4)oo’<2IB|a|-: 


n3/2slnh(^S) 

For  very  weak  fields  such  that  &L  «1>  we  have  flnally 

'hr ur 


(7.13) 


12YT 


- 1) 

i*  slnh(^SB) 


(7.14) 


The  last  term  in  (7.14)  arises  from  the  oscillatory  terms  of  (7.10) 
in  the  limit  as  3 H-^0,  We  see  that  it  is  much  larger  than  the 
true  steady  susceptibility  for  such  weak  fields  and 

sufficiently  low  t«mnp  ro  fnwo  W»s  + 4-V-t  ~ 4-  _ J • . 

• * w wAAcmu  xuo  0x5x1  axiu  iiidgai uuae 
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depend  critically  on  the  ratio  £ -fro. 

9 

We  note  from  (7.10)  that  the  ordinary  Landau  susceptibility 
is  the  dominant  steady  term  for  all  but  exceedingly  weak  fields. 
This  contrasts  with  Dingle's  finding  that  electrons  in  a box 
show  a steady  susceptibility  considerably  larger  than  the  Lan- 
dau value,-  for. weak  fields  comparable  to  those  required  for  the 
validity  of  (7*13)*  This  contrasts  also  with  the  importance  of 
the  "surface  diamagnetism"  for  the  box  model  at  intermediate 
fields.  However,  the  oscillatory  terms  in  (7.13)  are  quite 
similar  to  the  oscillatory  terms  in  Dingle's  weak  field  calcu- 
lations  if  we  identify  the  "radius"  R of  our  system  by  setting 
% - 


in  which  we  have  used  (5»9),  have  written  <o  * n . and  have  from 
(3.14)  h(0)“  § -9-sinOcosQ*  Here  rf(Q)  ig  cos2&sin9  or  cos20  sin2©. 
If  tf(0)  * coseO,  as  in  the  R2  term  of  (5.4)  in  (5.3),  then  since 
= - 2cos2®,  we  see  from  (5.13)  that  the  integral  is  exactly 
integrable.  In  general  we  must  resort  to  an  approximate  evalu- 
ation asymptotic  in  w.  The  following  procedure  does  not  yield 
the  complete  asymptotic  expansion  of  (Al)  but  gives  enough  of  it 
for  our  present  purposes. 

The  integral  in  (Al)  is  of  the  form  approximated  by  Steele2^ 
with  van  der  Corput's  method  of  critical  points. However,  the 
ultimate  justification  of  this  procedure  in  this  problem  seems  to 
be  the  method  of  steepest  descent,  since  for  n > m-1,  an  integral 
of  the  form  j xne*f^xmdx  does  not  converge  if  the  path  of  in- 
tegration Is  *'°  the  real  axis.  Since  we  shall  have  to  avoid 
divergences  resulting  from  terms  in  (am)  , with  k>  x+1,  in  the 
asymptotic  expansion  of  the  inner  integral  in  (Al)  when  the  u 
integration  is  carried  out,  we  should  justify  our  procedure  care- 


fully and  to  this  end  can  most  conveniently  use  the  method  of 
steepest  descent  directly.  The  following  discussion  is  tailored 
to  the  particular  rf(0),  h(0)  in  our  problem. 

t 

We  now  regard  «f( 0)  and  h(Q)  as  analytic  functions  of  the 
complex  variable  0.  From  the  imaginary  part  of  h(0)  we  find 
that  for  cou  > 0 the  path  of  integration  of  the  inner  integral  in 
(Al)  should  be  deformed  into  the  "steepest11  path  y1+  Y2  sketched 
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in  Fig,  VI.  The  complete  asymptotic  development  of  the  inner 

integral  is  then  obtained  by  expanding  rf(©)  and  h(0)  about  the 

critical  points  e where  h(Q)  assumes 

its  greatest  value  (unity)  on  the  pathyv*y0: 

1 2. 


ioU 

, TT 


i*l,2 


>t) ^ ^-7  ^(k)(01)C0-01)k^  ^ 


hCm)(o.)(o-ejm 


(A2) 


ifnU 

TT 


r-1. ^ t 1 Pi  da’ 


where  h ra  (Oj)  is  the  first  non-zero  derivative  of  hC©^  (for  our 
problem  the  third).  Each  term  in  the  expansion  is  evaluated 
asymptotically  (such  that  (tou)n  times  the  error  in  each  term 
approaches  zero  as  (u)U)  — >■  od  for  all  n)  by  replacing  by  the 
rectilinear  path  \l±  tangent  to  at  ©i  and  taken  in  the  same 
sense*  The  general  result  from  terms  with  p = 0,1  is  that  given 
by  Steele  in  his  equations  (A9)  and  (A10),® 

Using  the  special  form  of  h(0),  rf(0)  in  our  problem  (Table  IV), 
the  inter-- ting  terms  of  (A2)  are  those  involving  rf(2)(© ,) , ) 

and  /2)(©i)h  ’'(O^,  at  both  + 


Since  h^^(+^)  = -4  and 


nil 


(®  -5) 


TTNn 


- S>3  ~ (4?) 3 


r(ar> 

(A3) 


(U2) 

for  oju  > 0,  (with  the  corresponding  equation  at  - ~) , we  have  for 
the  first  two  terms  at  + for  the  inner  integral  in  (Al)  with 

rf(0)  * cos2©  sin©  and  tou  > 0 


* cos^©  sin©  and  o»u  > 0 

^ -fe5/3  e^r(5/3) 


g)5/3e^p(5/3)]. 


(kl) 
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Por  oju  < 0 we  obtain  the  development  by  conjugation,  as  with  the 
asymptotic  formula  (5°  14).  A similar  expansion  holds  for 
tf(  6)  - cos1 2&sin20. 

This  expansion  of  the  inner  integral,  which  we  shall  call 
♦(cm),  in  (Ai)  is  asymptotic  in  uju,  but  since  the  integral  over  u 
extends  down  to  u = 0 we  must  make  certain  that  the  result  after 
the  u integration  is  asymptotic  in  Since  a term  in  (l/u)k  in 
(A4)^  with  k£  x * 1,  loads  to  a divergent  u integral,  we  see, 
as  remarked  above,  ohat  we  can  at  best  obtain  by  this  method  no 
more  than  the  first  few  terms  of  the  complete  asymptotic  develop- 
ment of  the  double  integral  (Al). 

-We  accordingly  assume  an  asymptotic  approximation  to  ijr(am)  of 
the  form 

5 (am)  , (AJ) 

(tou) 

where  k < x + 1 and  |a(aju)|  < c,  a constant,  for  all  oju,  such 
that  for  an  arbitrary  6>0,  there  exists  an  A such  that  whenever 
o)U  > A 

|t(uxi)  - i;(um)|  (wu)^  8 , (A6) 

(The  proofnls  readily  extended  to  a more  general  expansion  of 
^(tou)  = > v » Provi8ed  that  all  the  k^  < x + 1. ) Then 

if  u >e  * ^ and  <o>  A/e  . (A6)  will  be  satisfied.  More- 

over for  either  of  the  4(Q)  of  interest  to  us,  j^(oni)|<  tt  . 
Accordingly, 


1 

| j (*<cou)-  j?(cou))ux(l-u)ydu 


u> 


< cokj J*  ijr(a)u)ux(l-u)y  du 


CO 


I j*  ^(cou)ux(l-u)ydu  + <ok  ^ iKcou)-  ^(cou) 


ux(l-u)ydu 


- ^1 


r E - 
Adu 


u 


♦ C 


re 

r 


u‘  ‘du  +6 


'1 


ux_k  du 
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TT U)^  X+l 

- £ 

X*1 


(X-k+TJ 


+(7^+TT 


By  choosing 


(A?) 


e 


where  A « •|(jt+l-k)  > 0,  (A8) 


we  may  determine  A such  that  <V3,  where  V .-C  is  arbitrarily 

small.  We  may  then  choose  co  such  that  the  first  two  terms  In  (A 7), 
which  according  tp  (A8)  each  vary  as  a negative  power  of  to , are 
each  less  than?2/3,  and  simultaneously 


x+l-k-X. 

V ^ A<  (A9) 

which  suffices  to  make  u)>  A/fe  . Hence  for  sufficiently  large  u> 


(a) 


J ^ - ^(o>u))  ux(l-u)ydu 


<n 


so  that  substitution  of  (A4)  yields  the  asymptotic  expansion  of 
(Al>  through  terms  with  k < x ♦ 1,  but  we  have  no  assurance  that 
higher  terms  are  given  accurately,  even  if  they  do  not  diverge. 
Substitution  of  (A4)  yields  then  for  (Al) 

^ Vx-l*y(a°  4 i(aC)5/3  • ^ p(5/3)  V ^ (u>) 

frrl  (A10) 

4 it  Vx-l*y(0)  *t(aS)5/3  • P ( 5/3 ) V * (0)  , 

and  we  must  use  (5.14)  to  eliminate  terms  varying  as  (l/cu)k  with 

^ ^ x 4 *or  we  see  that  this  condition  actually  elimi- 

nates all  the  oscillatory  terms,  so  that  the  oscillatory  term  in 
(5.28)  may  oe  inaccurate.  We  know  at  least  that  it  is  no  larger 
than  the  order  of  magnitude  there  given,  so  that  Table  I assures 

in  tho * 1ft  ^rtnlnl v 

— - - 

The  dominant  ter.  1„(A10)  1,  the  third  one,  but  this  vanishes 
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Appendix  g. 

In  order  to  calculate  the  second  term  of  (5.25)  and  the 

correction  to  i we  muj  Include  the  next  term*  - £§5| — (r_r  )3 

7 2mc2r0  0 

of  the  expansion  (3.9).  in  th<-  phase  Integral  (3.10).  We 

then  change  variables  to  x ard  y of  (3.13)  and  set  u « (ro_r)  (£fi) 

Making  the  further  substitutions 


v = u [ ± + 


* 1 + 


2©fR 


we  obtain  for  the  phase  integral 


2>±  (i  - 


) dv  * (n  + £ ) rr 


(Bl) 


2 

where  y * and  (Bl)  is  accurate  through  terms  in  1/R.  Intro- 
ducing the  parameter  z * v^yi1,  which  has  the  range  -l<z  <le  we 
obtain  the  equations  replacing  (3.14)  by  integrating g 


[2n+l)  n 


cos  Z-Zll-Z ) 


1+1 


#(2n+l>Jvir(l  - z2)3/2 


(||)i  Bfcos^z  -z(l-z2)i)3> 


(B2) 


(2n+l)fyn  | Tjztcos^z-zd-z2)^)-  ^(l-z2)$j 


(cos^z-zd-z2)*)*  (||)^R  dos”1z-z(l-z2)^)3/2 

(B3) 

Substituting  from  (3.13)  into  the  expression  for  (P-N^  obtained 
from  (5.3)  and  (5*4),  we  have  to  evaluate  an  integral  of  the  form 


j x(z)  (E  - 3Hy(z))*dz. 


(B4) 


We  cannot  expand  each  factor  in  a Taylor  series  in  1/R„  since  the 
last  factor  has  a branch  point  when  E * 3Hy(z),  and  we  wish  to 
integrate  over  z.  Hence  we  change  variables  in  (B4)  from  z to  t, 
where  z(t)  is  defined  such  that  with  accuracy  through  terms  in 
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y(z(t))  *=  yo(t), 


(B5) 


where  yQ(z)  is  obtained  from  (B2)  by  dropping  the  1/H  term  and 
is  the  same  as  in  (3.14).  We  find  readily  that  to  this  accuracy 
(B4)  becomes 


Zl<"^ 


ay  (t) 


(E  - 0Hyo(t))*  dt.  (B6) 


We  obtain  x(t+q(t)/R)  from  (B3)  after  determining  q(t)  from  (B5) 
and  (B2),  and  the  integral  (B6)  is  then  readily  evaluated  with 
the  methods  of  Section  V and  Appendix  A.  We  find  for  the  r = 0 
correction  to  (F-Nlp^ 


(.29)  ^ 1 V 

" D2  (PH)2  J0 


e7/2  dE 
dE  " ^ J 


(B7) 


which  cancels  the  corresponding  terms  of  (F-N £)g  an<^  (F-Ni£)  + in 
(5°  22)  and  (5»24),  in  agreement  with  Osborne's  theorem.  For  the 
correction  to  be  added  to  the  steady  part  of  (F-N^Jg  we  obtain 


OH) 


<%* 


(|)5/3  £(5/3X1  -7273) 


SI  V7T  T(  17/6) 


'ST  E 


11/6 


dE  , (B8) 


which  contributes  already  included  in  (5»29) 
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Appendix  C. 

To  prove  (6.8)  we  use  Poisson's  sum  in  its  usual  form^'L 


^ ]*«.) 

s=l  n=-  QD  4) 


e2nins  ds^  (C1) 


We  now  let  f(s)  = g(s+a),  0 < a < 1,  and  follow  the  method 
used  by  Dingle  for  a * 1/2: 

CP  CD  CD 

g(s**a)  * g(a)  + ^ j gCs+cOe2"1118  ds 
s=o  n=-  a)  ''-a 


-£  / 


(C2) 


g(s+a)e2Tlins  ds. 


The  general  Fourier  series  for  a function  defined  in  (-L,L)  is 


od  - ninx 

p(x)  ~ ii  e L 

n=-  od 


L 


nln*' 


p(x')  e 


dx 


(C3) 


in  which  the  series  converges  to 
of  ordinary  discontinuity  and  to 
Hence  setting  L = ^ and  defining 

g(x  ♦ a - for  3 -a  < x < - 


p(x-O)  + p(x+0)  | at  a point 
p(L-O)  + p(-L+0)J  at  +L. 


p(x)  = 


2 2 ’ 
for  - i 1 


(C4) 


<x  - a, 


we  find  on  substituting  in  (C3)  and  setting  x = 


= 1 


_ge_ 


£«<«>.  2_vi r •i)«2nin5 1* 


n=-  oo 
. 00  /p 

J 

n=-  od  J- a 


•t- a 
g(s+a)e2nins 


(C5) 


ds, 
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Hence  the  first  and  third  terns  on  the  right  of  (C2)  cancel,  and 
we  have  left 


g(s*a)  * e-2"1"0  f g(s)  e2"1*3  ds.  (C6) 

s=o  n=-  cd  Jo 
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17.  Since  V(r)  has  an  infinite  step  at  r = R,  we  should  impose  on 
f(r)  the  boundary  condition  frR;  = 0.  This  leads  in  (3*5) 

to  (n  + 3/4)  in  place  of  (n  + 1/2)  when  the  upper  limit  of 
integration  is  R,  although  it  is  not  clear  at  what  point  as 
r2  approaches  R we  should  start  to  vary  from  1/2  to  3/4.  This 
replacement  does  not  alter  the  results  of  Section  TV  or  the 
derivation  of  Mo  in  Section  V but  does  change  the  correction 
terms.  This  matter  is  discussed  at  length  in  Section  VI,  but 
until  that  section  we  shall  use  (n  + 1/2).  ?rope~  use  of 
(n  + 3/4)  would  cause  the  curves  of  Fig.  I to  rise  at  a 
slightly  lower  value  of  x and  remain  slightly  above  the 
curves  plotted.  It  would  undoubtedly  smooth  out  the  low 
end  of  the  moment  curves  in  Fig.  II. 

18.  Peirce  . A Short  Table  of  Integrals.  New  York,  1910,  No.  187. 

19.  The  effect  of  spin  should  be  determined  b Jr  Dingle’s  procedure 
(see  reference  3,  I).  This  adds  to  the  steady  moment  the 
usual  spin  paramagnetism  and  changes  the  phase  of  the 
oscillatory  terms* 

20 o The  tail  is  clearly  defined  only  if  we  use  approximate 

methods,  such  as  the  WKB  procedure,  to  determine  the  energy 
eigenvalues.  In  an  exact  solution  all  eigenvalues  would 
depend  on  ro,  hut  the  departure  from  the  free  electron  eigen- 
values would  Ke  appreciable  only  for  states  with  ro  ^ R. 

21.  When  we  obtain  M from  (5«3)we  should  have  no  terms  Involving 
[38(n.H,E)]/dH  or  dAr  _/dK  and  in  fact  do, not,  for  although 
both  8(n,H,Ej  and  Ar0n<iepend  on  H,  (E-e)t  vanishes  at  r =6 
and  de/dro  vanishes  at  ro=R  - Arori9  so  that  the  integrand  is 
zero  at  each  limit  of  integration. 

22.  E.  T.  Whittaker  and  G.  N.  Watson,  A Course  of  Modern  Analysis, 

4th  ed.  Cambridge,  1952,  p.  337 . 

23.  Whittaker  & Watson,  £E>.  ai±.  Mvm(z)  is  incorrectly  given 
in  terms  of  Wk  m(z)  and  W-k  m(“z)  in  Ex.  2 offci6.4l,  p.  346. 

The  expression’ should  be  corrected  by  replacing  ekm.  by  e-krri 
in  each  term.  The  asymptotic  exDansion  of  Wk  m(z)  is  given  on 
P.  343. 

24.  This  theorem,  one  form  of  which  was  previously  suggested  by 
Dingle  (Ref.  3,  III),  has  proved  to  be  true  to  the  accuracy 
used  In  every  problem  in  which  it  has  been  tested,  including 
our  exact  solution  of  the  Darwin  problem  (Section  VII).  Ko 
one  has  as  yet  given  a general  proof  $ Osborne's  proof  is 
valid  for  fairly  general  potentials  but  requires  the  use  of 
the  WKB  approximation,  which  is  known  to  he  accurate  only 
asymptotically  as  either  (fi2/m) — *~0,  or  as  the  quantum 

number  n becomes  large  (J.  1.  Dunham,  Phy*.  Rev.  41,  713  (1932)). 
Since  this  is  in  this  sense  a semicl<assical  approximation,  it 
is  not  impossible  that  its  use  obscures  a purely  quantum  effect 
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not  shown  by  Darwin’s  model. 


25.  The  magnitude  of  the  Mn  corraction  was  incorrectly  given  in 
the  abstract  of  an  earlier  report  on  this  work.  Moreover, 
at  the  time  of  that  report  the  work  of  Section  YI'had  not 
been  done.  (Bull,  Amer.  Phvs.  Soc.  28,  no.  1 WA8.) 


2fc.  Several  of  our  small  oscillatory  terms  correspond  to  terms 
found  by  Osborne.  Thus  the  second  term  in  his  M-i  0 0 is 
of- the  same  form  as  the  oscillatory  part  of  our  M{,’and  his 
■i  1 © Is  like  the  oscillatory  part  of  M+  , or  M2  and  arises 
frdm’the  energy  surface  near  states  with  positive  s, 

2 7.  L.  I.  Schiff,  Quantum  Mechanics . McGraw-Hill.  New  York, 
1949,  p.  186. 

28.  Whittaker  and  Watson,  £2.  cit. , p,  244. 


29.  See  Reference  3,  Appendix.  Steele’s  integrals  are  very  simi- 

lar tc  ours.  Thus  his  In  is  apart  from  a constant  the 
imaginary  part  of  our  Vo  -1  (u>/2).  His  Ii2  is  proportional 
to  Im  Vx}-^.(w/2),  and  hil  ?I2i+I22  rl 

* 1 f dz 

^ 1 

30.  Proceedings  of  the  Section  of  Sciences  T Nederlandsche  Akademie 
Van  Wetenschappen  5l,  650  (194-8)7 


31o  Reference  3,  I (Appendix).  The  conditions  for  the  validity 
of  the  formula  and  the  sense  in  which  the  convergence  of 
the  series  is  to  be  understood  are  obtained  from  Titchmarsh, 
Introduction  to  the  Theory  of  Fourier  Integrals , Oxford, 
I932,  p.  60. 


